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Abstract 

A variational principle is introduced which minimizes an action formulated for con- 
figurations of vacuum Dirac seas. The action is analyzed in position and momentum 
space. We relate the corresponding Euler-Lagrange equations to the notion of state 
stability. Examples of numerical minimizers are constructed and discussed. 
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1 Introduction 



In the standard model of elementary particle physics many structures and parameters 
are built in ad-hoc. In particular, the theory makes no statement on the masses of the 
leptons and quarks; the mass parameters need to be put in empirically from experimental 
data. This situation is not quite satisfying, and one would hope for a more fundamental 
explanation for the ratios of these masses. As an approach towards a theory which makes 
predictions for the mass parameters, we here set up a variational principle where we 
minimize an action formulated for configurations of vacuum Dirac seas. 

The minimization problem is nonlocal, nonlinear and involves two or four free param- 
eters. Therefore, a systematic study of the structure of all minimizers goes beyond the 
scope of this paper. But we construct and discuss numerical examples of minimizers. The 
technical core of the paper is to compute the Fourier transform of the action and to relate 
the regularization on the light cone to a suitable regularization procedure in momentum 
space. This transformation to momentum space simplifies the numerical analysis. Fur- 
thermore, it reveals a connection between the corresponding Euler-Lagrange equations and 
the notion of state stability as introduced in [TJ §5.6] for vacuum Dirac sea configurations. 

Since the action principle for vacuum Dirac seas seems of interest by itself, this paper 
is self-contained and easily accessible. Nevertheless, our action principle can be seen 
in a broader context as a modification of the variational principle introduced in [1] for 
the fermionic projector in discrete space-time (see also [2] or the review articles O d]). 
More specifically, we modify the action of pQ such as to make it finite in the Lorentz 
invariant setting, treating the singularities on the light cone by a suitable regularization 
procedure. Our method is complementary to that in [H Chapter 4]: whereas in pQ the 
behavior of the singularities on the light cone is considered in the so-called continuum 
limit, we here disregard these singularities completely and concentrate instead on the 
regular contribution away from the light cone. Yet another method is used in [5], where 
the vacuum Dirac seas are regularized, and it is analyzed how the operator product M.P 
behaves as the regularization is removed. These different approaches all complement each 
other and should eventually be combined to get a full understanding of the structure of 
the minimizers of our variational principle. For an overview of the different methods we 
refer to the review paper [6]. 

The paper is organized as follows. In Section [2] we set up our variational principle 
in position space and discuss the relation to the variational principle in discrete space- 
time as introduced in [I]. In Section [3] we prove Plancherel-type formulas for Lorentz 
invariant functions. These Plancherel formulas are basic for the subsequent reformulation 
of our variational principle in momentum space. In Section [4] we outline the basic method 
and explain the connection to state stability, whereas Sections [5] and [6] are devoted to 
the more technical and computational aspects. In Section [7] we introduce two additional 
free parameters, which arise in the context of state stability, and built them into our 
framework by introducing the so-called extended action. In Section [8] we outline our 
numerical methods and discuss a few examples of numerical minimizers and state stable 
Dirac sea configurations. 

2 A Variational Principle for Dirac Sea Configurations 

In this section we shall formulate a variational principle for Lorentz invariant Dirac sea 
configurations. For clarity and self-consistency, we proceed without referring to the vari- 
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ational principle in discrete space-time as introduced in [T]. The connection to the latter 
variational principle will be explained later in this section (see after Definition 12 .2|) . We 
begin with a general ansatz for the fermionic projector of the vacuum involving an ar- 
bitrary number of generations g G N of Dirac seas of masses mi, . . . ,m g G R, which are 
taken into account with weight factors pi, . . . ,p g > (for a discussion of the weight factors 
see [5j Appendix A]). More precisely, we introduce the distribution in momentum space 

9 

P(k) =Y,Pp(Jf + m p ) S(k 2 - m}) Q(-k°) (2.1) 
13=1 

and consider its Fourier transform to position space 

P(0 = / (04 P (k) . (2.2) 

Here k£ denotes the inner product of signature (H ) in Minkowski space M. 

We first define our action formally and give the analytic justification afterwards. For 
the upper and lower light cone and its boundary we use the notations 

/ v = {£ \e > o, e° > o} , / a = {£ i e > o, e < 0} , l = ku 2 = o}. 

For any £ G I v we introduce the so-called closed chain A and its trace-free part Aq by 

A = P(C) P(0* . A) = ^ - i Tr(A) , (2.3) 

where the star denotes the adjoint with respect to the indefinite inner product on the 
Dirac spinors (where ^ = is the adjoint spinor). Thus 

where the dagger denotes the transposed complex conjugate matrix. Due to Lorentz 
invariance, Aq can be written in the upper light cone as 

M = \ f{f) for Z G / v , (2.4) 
and the relation ^4q = implies that the function / is real. Hence the Lagrangian 

C = Tr(Al) = en?? (2.5) 
is non- negative and depends only on the positive parameter z := £ 2 . Thus the integral 

formally , ^ % ^ ^ 

formally defines a positive functional depending on the free parameters mi, . . . ,m g and 

Pl,---,Pg- 

Let us now give the definitions (|2.3l 12.51 12.6P a rigorous mathematical meaning. An 
explicit calculation of the Fourier transform (|2.2p (see [5j Section 3]) shows that the dis- 
tribution P(£) has the following properties. First, it is singular on the light cone L, but 
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is a smooth function otherwise. Furthermore, it is Lorentz invariant. This means in the 
upper light cone that P(£) can be written as 

Pit) = iv(?) + h(?) for£e/ v (2.7) 

with smooth complex functions v,h £ C°°(]R + ), which can be given explicitly in terms of 
Bessel functions. Thus the pointwise computations (|2.3p and (|2.5[) make sense and 



f(z) = Re(v(z)h(z)) G C°°(l 



We conclude that the integrand in (|2.6p is a smooth function. Near infinity, an asymp- 
totic expansion of the Bessel functions yields that the function / decays like 0(z -2 ), and 
thus the integral in (12.6f) is absolutely convergent at infinity. Near the origin, / has the 
expansion (see O eq. (3.5)]) 

f(z) = ^ + ^ + O(logz), (2.8) 

z z z 

where we set 

9 



1 

m 3 = Yl PaPp( m l + m l) ( 2 - 9 ) 



a,/3=l 
9 



1 

™ 5 = no 5 y2 Pa P(3 (m a - nifs) 2 (m a + nip) 3 . (2.10) 

a, /3=1 

Thus according to (|2.5p . the integrand in ()2.6[) has a non-integrable pole at z = 0, 

Tr(^) z = /(z) 2 z 2 = ^ + + ©(log z) . (2.11) 

z z z 

In order to make sense of the integral, we need to subtract suitable counter terms. The 
simplest method is to add indefinite integrals of the poles in (I2.1ip evaluated at z = e, 

lim (J Tr(A§) z dz - ^ + 2 m 3 m 5 loge^j . 

However, since the indefinite integrals are determined only up to constants, it seems natu- 
ral to allow for adding arbitrary constants times rri3 and msms. This freedom of modifying 
the action will turn out to be very useful later. In order to be even more flexible, it is 
preferable to treat the freedom more generally by a function i 7 (1113,1115). This gives rise to 
the following definition. 

Definition 2.1 (Lorentz invariant action) For any given function F G C°°(IR x R, R) 

we define the action S = 5 (mi, . . . , m g ,pi, . . . , p g ) by 

( f°° m 2 \ 

S = Jimfy Tt(A 2 ) zdz - -± + 2m 3 m 5 log e\ + F(m 3 ,m 5 ) . (2.12) 



Here Aq is defined by &2.3\) for any £ G I v and z = £ 2 . The parameters ITI3 and ms are 
defined by (KM \KM- 
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Note that due to the negative counter terms, the action (|2.12|) need no longer be positive. 
But for any fixed e > it is clearly bounded from below. We remark that in Definition 17. II 
we shall extend this action by additional summands involving two more free parameters. 

Next we want to introduce a corresponding variational principle. In order to avoid 
the trivial minimizers p\ = ■ ■ ■ = p g = or mi = ■ ■ ■ = m g = 0, we need to impose a 
constraint. 

Definition 2.2 (Lorentz invariant variational principle) We minimize the action S, 
112. 12\) . varying the parameters pi, . . . , p g > and mi, . . . , m g > under the constraint 

g 

T:=Y,Ppr4 = I. (2.13) 

p=i 

In order to understand the relation between the above variational principle and the 
variational principle introduced in discrete space-time in [IJ §3.5] and [2], we first consider 
the Lagrangian (|2.5p for vectors ^ which are space- like or inside the past light cone I A . 
Taking the adjoint of (|2.2p and (|2.ip . one sees that 

p(-0 = p(0* ■ 

Hence the closed chain (12.3|) can be written as A = P(0 P(— £). Due to the ansatz (12. H 
12. 2p . the matrices -P(±£) are of the form a±$ + (3± with ct±,f3± G C. This shows in 
particular that the matrices -P(£) and P(— £) commute, and thus 

A(-0 = A(C) and A (~0 = MO ■ ( 2 - 14 ) 

For space-like ^, we can use Lorentz symmetry to write the matrix -Ao(£) in the form 
MO = (d{i 2 ) with a function g G C°°((-oo, 0)). Replacing ^ by and using ([2TT4]) , 
we see that g must vanish identically. Combining this with (12. 41) and (|2.14p . we conclude 
that away from the light cone, Aq is of the form 

MO = |/(e 2 )e(C 2 )e(e°) if^L (2.15) 

(where e is the step function e(x) = 1 if x > and e(x) = — 1 otherwise). In words, 
is reflection symmetric (|2.14|) . causal and Lorentz invariant. Moreover, for timelike ^ 
the roots Ai, . . . , A4 of the characteristic polynomial of A (counted with multiplicities) are 
computed to be real. According to [5j Lemma 2.1], these roots all have the same sign. 
Combining these facts, we can write the the Lagrangian (12.51) in the alternative form 

C = Tt(A 2 ) - i Tr(A) 2 = \A 2 \-^\A\ 2 if £V 0, 

where | . | is the spectral weight (see [TJ §3.5]) 

\A\ = ^N. 

i=l 

This is precisely the Lagrangian in the "model example" introduced in [TJ §3.5], for a 
particular value of the Lagrange multiplier. The same Lagrangian also appears in [2] 
as the so-called critical case of the auxiliary variational principle. As this Lagrangian 
vanishes for space- like £ and is symmetric under the transformation £ — ► — £, it was no 
loss of generality in Definition 12. II to restrict attention to £ in the upper light cone. Hence 
the only differences between the above variational principle and the variational principle 
in [H [2] are the following points: 
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(i) A regularization on the light cone is used in (|2.12p . 

(ii) The additional constraint (|2.13j) appears. 

(iii) The method of summing or integrating over space-time is different. We can describe 
the transition from the action in [TJ [2] to (|2.6|) by the replacement 



E • 

x,y£M 



z dz 



Let us briefly discuss these points. It is not astounding that a regularization on the light 
cone is necessary, because the fermionic projector has poles and singularities on the light 
cone. Such regularizations have been used previously in the derivation of the continuum 
limit [U Chapter 4] and for the vacuum fermionic projector in [5j. The new and pleas- 
ant feature of the regularization (i) is that it does not break Lorentz invariance. The 
necessity of the constraint (ii) can be understood as replacing the condition in discrete 
space-time that the number of particles is fixed. However, the detailed form of the con- 
straint (I2.13|) will not become clear before Section HI The replacement (iii) requires a 
detailed explanation. The natural method to describe the transition from discrete to a 
continuum space-time is to replace sums by space-time integrals, 



E 

x,y&M 



d 4 x / d 4 y- = / d 4 x / A' 

M JM JM JM 



where in the last step we made a change of variables. For homogeneous systems as consid- 
ered here, where the fermionic projector depends only on the difference y — x, the integrand 
of the rr-integral is a constant. Thus this integral merely gives an infinite constant, and it 
seems natural to simply drop it. Thus it remains to discuss the replacement 

/•oo 

■•• A — ► / ■■■zdz. (2.16) 

M JO 

This replacement has no simple justification; it changes the variational principle substan- 
tially. The integral on the right has the nice property that it can be finite for Lorentz 
invariant integrands, whereas the left integral is necessarily infinite for Lorentz invariant 
integrands, because the integrals over the hyperbola £ 2 = const diverge. Furthermore, 
the integration measure z dz in (|2.16p can be motivated from the fact that it has the 
same dimension of length four as the measure d 4 £. However, we again point out that 
the replacement (|2.16p changes the variational principle completely. Thus the variational 
principle introduced here is not merely a "continuum version" of the variational principle 
in [HE]. It is a different variational principle. Nevertheless, as we shall see in Section [H 
there is a close connection between minimizers of the variational principle of Definition ^. 2| 
and the notion of state stability as introduced in [TJ §5.6]. 

We close this section with the derivation of the corresponding Euler-Lagrange (EL) 
equations. We denote the first variation of the masses mp and weights pp by 5mp and 5pp, 
respectively. We let 5P be the corresponding first variation of the fermionic projector (|2.1l 
12 . 2|> . The variation of the Lagrangian is 

5Ti{Al) = 2Tt{A q 8Aq) = 2Tr(A 8 A) , 
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where in the last step we used that the scalar component of SA drops out of the trace 
because the other factor Aq only has a vector component. Using f|2.3|) and (|2.14p . we 
conclude that 

6Tr(A%) = Tr(A (5P(t)P(0* + P(0SP(tT)) = 2ReTr(A P(£) SP(~o) ■ 
Treating the constraint with a Lagrange multiplier A, the EL equations become 

lim (J °° Re Tr(^ P(0 SP(-£)) zdz - + {5m 3 m 5 + m 3 <5m 5 ) log ej 

+5F(m 3 ,m 5 ) - X5T = 0. (2.17) 

This equation looks rather complicated, and in this form it does not seem helpful for 
understanding our variational principle. In order to clarify the situation, in the next 
sections we shall transform our variational principle to momentum space. 



3 A Plancherel Formula for Lorentz Invariant Causal Func- 
tions 

Studying the Lorentz invariant action (|2.12j) directly in position space has the following 
disadvantages. First, expressing Aq explicitly using Bessel function, one sees that Aq(z) is 
oscillatory for large z, making it difficult to analytically estimate or to numerically compute 
the integral in (j2. 12|) . Second, the corresponding EL equations in position space (|2.17p have 
no obvious interpretation. For these reasons, it is preferable to transform the action (|2.12p 
to momentum space. To this end, we need a Plancherel- type formula, which relates the 
integral over the position variable z = £ 2 to an integral over a momentum variable a := k 2 . 
In this section we shall derive such Plancherel formulas. For clarity, we begin with scalar 
functions. Thus, leaving out the Dirac matrices in (|2.15p . we consider functions of the 
form 

F(0 = /(ae(£ 2 )e(£°), G(0 = <7(£ 2 )e(£ 2 )e(£°) (3.1) 

with measurable complex- valued functions / and g. In view of the integration measure 
in (12. 6p . we introduce on such functions the inner product 

(F,G) := / f(z)g(z)zdz. (3.2) 
J 

We denote the space of functions for which the last integral converges absolutely, by 
L 2 (M,zdz). For the upper and lower mass cone and their union we introduce the notations 

C v = {k\k 2 > 0, /c° > 0} , C A = {k I k 2 > 0, k° < 0} , C = {k \ k 2 > 0} . 

Taking the usual Fourier transform in Minkowski space, 

F(k) = J dV(()e^, (3.3) 

a symmetry argument shows that F(—k) = —F(k), and from Lorentz invariance we 
conclude that supp-F C C. Hence F can be written similar to (|3.ip as 

F(k) = f(k 2 )0(k 2 )e(k°). (3.4) 
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On such functions we introduce similar to (|3.2p the inner product 

(F,G) := ^ J™J(a)g(a)ada, 

where we set a = k 2 . The space of functions, for which the last integral converges abso- 
lutely, is denoted by L 2 (M,ada). 

Theorem 3.1 (Lorentz invariant Plancherel formula, scalar case) For functions 
of the form $3. the Fourier transform 113. 3\) is a unitary mapping from L 2 (M, zdz) to 
L 2 (M,ada). In particular, for all F, G G L 2 (M,zdz), 

1 



f(z)g{z)zdz = -j—^ J f(a)g(a)ada. 

Proof. A possible method would be to express the transformation from f to f with Bessel 
functions, and to use properties of the so-called Hankel transformation (see for example [10] 
and |9j). In order to avoid working with special functions, we shall here use a different 
method, which also has the advantage of explaining the integration measures z dz and a da. 
Nevertheless, we briefly outline the connection to the Hankel transformation: Rewriting 
the Fourier integral (|3.3p as 

F(k) = / f(z) dz / A 5(e - z) e(£°) e"** , 



we can carry out the last integral using Bessel functions. Comparing with (|3.4|) we obtain 

f(a) = 2nr / f(z) a dz . 

Jo V az 

After transforming to the variables x = \[z and y = -^/a, we can apply Parseval's equation 
for the Hankel transform [10\ Theorem 5.1.2]. 

We now give the alternative proof which avoids special functions. Using a standard 
approximation argument, it suffices to prove the theorem for functions F and G of the 
form (|3.ip with smooth and compactly supported f,g € Cg°(M+). A short calculation 
shows that then the wave operator becomes 

-□F(f) = (Wm 2 ) e(£°) , 

where W is the ordinary differential operator 

4 d ( 2 d 



Integrating by parts, one sees that W is a non-negative symmetric operator on the Hilbert 



space L 2 (M. + , zdz) with domain of definition 

In order to construct a self-adjoint extension of W, we use the well-known fact that 
the operator —d 2 is self-adjoint on L 2 (M. + ,dz) with domain V(—d 2 ) = Hq' 2 (M. + ). We set 

V(W) = {u measurable with zu(z) & T>(—d 2 )} . (3.6) 

To verify that W is self-adjoint, we use the definition of self-adjointness. Suppose that for 
given v,w G L 2 (R + ,zdz), 

(Wu,v) L 2 {R+tZdz) = (u,w) L 2 {R+jZdz) WueV(W). (3.7) 
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A short calculation shows that 

Wu = —d 2 u with := zu(z) . 
This allows us to rewrite (|3.7p as 

(-d 2 z u,v) L 2 {R + :dz) = (u,w) L 2 {R + 4z) \/ueV(-d 2 z ). 

Since the operator (—d 2 ) is self- adjoint, it follows that v £ T>(—d 2 ) and that —d 2 v = w. 
Hence v 6 T>(W) and Wv = w, showing that W with the domain (|3.6p is indeed self- 
adjoint. 

The spectral theorem yields 

(f,9)L 2 (9+,zdz) = / (/> di£aff)z,2( R + )/Zd2 ) ) (3.8) 

where cZ£/ a is the corresponding spectral measure. In order to get more information on 
the spectrum, it is most convenient to work in momentum space, where the operator — □, 
and consequently also the operator W, coincide with the multiplication operator f(k) i— > 
k 2 f(k). Then it becomes obvious that a(W) = K + U {0} and that the spectral measure 
of W is absolutely continuous with respect to the Lebesgue measure da. Furthermore, the 
functional calculus can be expressed by 

(HW)f)(b) = h(b)f(b), 

implying that the spectral measure satisfies the relation 

(dEJ)(b) = 5{b-a)f{b)da. (3.9) 
As a consequence, the integrand in (|3.8j) can be written as 

(f,dE a g) L2{R+:Zdz) = f{a)g{a) p(a) da (3.10) 

with a non-negative measurable function p. 

In order to determine p, we compute the left side of f|3. 10|) with the help of (|3.9p to 
obtain 

f'°° f d^k 

f(a)g(a)p(a) = j z dz — ^ f(h?) @(k 2 ) e(k°) 

X / Pue^ 9 (k> 2 )6(k' 2 -a)e(k>°), 
Jm ( 27r ) 

where uj and uj' denote the zero components of k and k', respectively. Using the scaling a — > 
X 2 a, k — ► Xk, k' — > Xk' and z — > z/X 2 , we find that 

f(X 2 a)g(X 2 a) p(X 2 a) = f(X 2 a)g(X 2 a) X 2 p(a) 

and thus p(a) = ca with a positive constant c. Using this fact in (|3.10p . the abstract 
formula (|3.8p becomes 

poo 1*00 

/ f(z)g(z)zdz = c f{a)g(a)ada. 
Jo Jo 
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We finally determine c using the symmetry between position and momentum space to- 
gether with the fact that the definitions of the Fourier transformation and its inverse 
differ by factors of (2-7r) 4 (compare (|2.2p and (|3.3p ). ■ 



We now extend the last theorem to include Dirac matrices. We thus consider similar 
to (|3.ip a function F (and analogously G) of the form 

f(0 = f /(£ 2 ) ®(£ 2 M£°) ■ (3.11) 

In the inner product we combine the two resulting factors £ to a factor z, 

(F,G) : = / zf(z)g(z)zdz. (3.12) 
J o 

Then the Fourier transform F(k), again defined by (|3.3H . can be written as 

F(k) = lf(k 2 )@(k 2 )e(k°). (3.13) 

The inner product in momentum space is defined by 

„ „ i r°° - — 

{F,G) = J af(a)g(a)ada. (3.14) 

With a slight abuse of notation, we denote the Hilbert spaces corresponding to the inner 
products (|3.12p and (I3.14|) again by L 2 (M,zdz) and L 2 (M,ada), respectively. 

Corollary 3.2 (Lorentz invariant Plancherel formula, vector case) For functions 
of the form $3.11\) . the Fourier transform \3. 3\) is a unitary mapping from L 2 (M, z dz) to 
L 2 (M, a da). In particular, for all F, G G L 2 (M, z dz), 



zf(z)g{z)zdz = -j—r i J af(a)g(a)ada. 
Proof. We can again assume that f,g£ C^°(M+). Introducing the scalar function 

f s (o = f s (e)®(n 2 )e(e) with /„=/, 

we obviously have F(£) = $F s (£)/2. Since multiplication in position space corresponds to 
differentiation in momentum space, we know that 

F(k) = l -^F s (k). 

Representing F s in the form (|3.4p with an overall subscript s, it follows that the function / 
in (13.13P is given by 

f(a) = 2f' s (a). 

Finally, we know that £ 2 -F s (£) = — □fc-Fs(fc) and thus 

zj^z) = Wf s (a) , 
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where W is minus the wave operator in momentum space, i.e. in analogy to (|3"T 



4 d ( 2 d 



W-X')--;sKi''l' (3 ' 15) 



Using for G the same notation as for F, we obtain 



z f(z)g(z) zdz = / z f s (z) g s (z) z dz 
Jo 

(*) 1 f°°^m ~ / s , 1 



zf s {a)g s {a)ada = — — j / ^^(^^^(a) atk 



, (27r)4 y o 

(+) 1 f°„!^T,/, s 2 , 1 



(2tt) 4 7 4 ^( a )5s( a ) a da = T^jAi y o af(a)g(a)ada, 

where in (*) we applied Theorem 13. 1\ and in (+) we used the explicit form of W, (|3.15|) . 
and integrated by parts. ■ 



4 The Variational Principle in Momentum Space, 
Connection to State Stability 

Using the Lorentz invariant Plancherel formula of Corollarv l3.2l we shall now transform our 
action and the corresponding EL equations to momentum space. The main difficulty is to 
treat the singularities of Aq on the light cone. We here merely describe the regularization 
schematically; the detailed constructions involving a rather subtle regularization technique 
will be carried out in the subsequent Sections [5] and [6l 

We want to apply Corollary !3.2l to the integral in (|2.12p . The simplest approach would 
be to regularize Aq by multiplication with a characteristic function, 

AKO ■= e(e-e)A (Z). (4.1) 

According to <|2.15j> . 



AUO = §/ e (*WM*°) with f{z) = ®{z-e)f(z) 



Then 



/OO /'OO 
Tr (A%) zdz = J Tr((^) 2 ) 



z dz 



f(zfz 2 dz = -j—^J \fc(a)\ 2 a 2 da, (4.2) 

where in the last step we applied Corollary 13.21 Here f e is defined via the relation 

J d\ 9(f - e) A = | f e {k 2 ) &(k 2 ) e(k°) . (4.3) 



The problem is that in (|4.3|) we take the Fourier transform of a function which is dis- 
continuous on the hypersurface {£ | £ 2 = s}. As a consequence, the function f £ depends 
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on the regularization in a complicated oscillatory way, making it difficult to analyze the 
e-dependence of the last integral in (|4.2p . 

To bypass this difficulty, we will use a different regularization method. First, we ex- 
tend Aq in the distributional sense across the light cone. More precisely, we introduce M (£) 
as in [5] by the requirements that 

M£S'(M) and M(0 = 2A Q (£) V£ L , (4.4) 

where S'(M) denotes the space of tempered distributions. Furthermore, M. should respect 
the symmetry (|2.14p . 

M(-£) = M(0- (4-5) 

The Fourier transform of M is denoted by M. G S'(M). Using again the symmetry 
argument after (|3.3p . we find that 

suppMcC and M(-k) = M(k) . (4.6) 

Next, we regularize A4 in a Lorentz invariant way which is simple and explicit both in 
position and momentum space. The regularized M is in L 2 (M,zdz), making it possible 
to apply Corollary 13. 21 This allows us to rewrite the action (|2.12p as 

S = ^ ^(J^Tr(M e (k) 2 )ada + F £ (m 3 ,m 5 )) , (4.7) 

where k £ M is any vector with k 2 = a. Here the new function F £ also involves the counter 
terms. For the first variation we then obtain 

5S = (2^4 S3 Of \ Tr s ^ £ ) ada + 6F ^ m ^) ■ ( 4 - 8 ) 

According to (|4.4p and (|2.3p . the distribution _M(£)/2 coincides away from the light cone 
with the trace-free part of the product -P(£) times P(— £). Since multiplication in position 
space corresponds to convolution in momentum space, one might expect that A4 equals 
the convolution of P by itself. More precisely, the change of variables 

»(p-9)£ 



suggests that it should be helpful to decompose M. £ as 

7W^(A;) = JF(/c) + AA e (fc) , (4.9) 

where we set 

/j4 
— 4j (P(ifc + g) P(g)") , defined if k 2 > 0. (4.10) 
(27TJ 4 V /0 

Here the subscript zero means that we take only the vector component, and the new 
function N £ involves the dependence on the regularization. As the poles of M are char- 
acterized by rri3 and ms (see fj2.8|> ) . it is not astonishing that J\f e will depend only on the 
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Figure 1: The convolution integral in (|4.10p . 



parameters rri3, rns and e. Using that P is supported on the lower mass shells (see f|2. 1[) ) . 
one sees that the integrand in (|4.10p is compactly supported (see Figure [TJ in the case 
when k £ C v ; the case k £ C A is similar). As a consequence, the integral in (|4.10p will be 
finite. Since M £ — ► M as e \ 0, it is obvious from (|4.9p that N £ converges in this limit 
to a function N(k) defined if k 2 > 0. Using (|4.9p in (|4.8p . we obtain 

5S = — \j lim / ada ^ Tr (M £ {k)5F{k)) (4.11) 



(2vr) 4 £ \o7 2 
+^ lim QH 1 Tr (.M £ MA e ) a da + 5F,(m3,m 5 )) . (4.12) 

We shall see that in (|4,lip the regularization can be removed. In (|4.12p . on the other 
hand, the integral diverges as e \ 0, but this divergence is compensated by the counter 
terms in 5F e . We thus obtain the following result, where for notational convenience we 
omit the irrelevant prefactors (2ir)~ 4 . 

Lemma 4.1 The first variation of the action h2.12\) can be written as 

SS = a da J —^Tr(M(k)6(P(k + q)P(q))^ + oF(m 3 ,m 5 ), (4.13) 

where k £ M is any vector with k 2 = a, and 5F denotes the variation of the function F 
in Definition \2.1[ 



This concludes our sketch of the regularization procedure. This sketch will be our guideline 
for the constructions in Sections [5] and [6l The proof of Lemma 14.11 will be given at the 
end of Section El 

We next consider the variation of the product in (|4.13p . Changing variables 



/ 



d A q (6P{k + q)) P(q) = J d 4 q (5P(q)) P(-k + 



and using that commuting the two factors in the integrand does not change their vector 
component, we find that varying the first factor is the same as varying the second factor 
if we flip the sign of k. In view of (14. 6h we can thus write the variation of the action as 



5S = J ada J -j*TrQM(k)P(k + q) + M(-k) P(-k + qj\ 5P{qi 

+5F(m 3 ,m 5 ), (4.14) 



13 



Figure 2: The g-integral in (|4.15p in the cases k € C v (left) and k £ C A (right). 




-k° 



-k 



supp P(k + q) 




k = 



Figure 3: The corresponding fc-integrals in the two cases of Figure [2j 



where k G M is an arbitrary vector with k 2 = a. Note that we here omitted the subscript 
zero. This does not change the value of the trace, because the factor A4 only has a vector 
component. 

The following geometric argument is crucial for getting an understanding of the EL 
equations. It suffices to consider the first summand in the square bracket in (|4.14p . because 
the second term is then obtained by flipping the sign of k. Using the support property 
in (14. 6h . we may insert the integral over a ^-distribution, 



a da / d 4 qM(k)P{k + q)5P{q) 



oo poo 

a da db I d*q 5(q z - b) M (k) P(k + q) 5P(q) . 
o Jo 



(4.15) 



For fixed a, b and k S C, the (/-integral extends over a compact subset of the hyper- 
bola {q\q 2 = b} (see Figure [2]). Any two vectors on this hyperbola can be transformed 
into each other by a Lorentz boost. Thus, using Lorentz symmetry, we can assume that 
the vector q is fixed. Instead, the vector k will no longer be constant, but it will run over 
a compact subset of the hyperbola {k \ k 2 = a} (see Figure [3]). Taking into account the 
relative sizes of the g-subset and the corresponding /c-subset, we obtain the transformation 
rule 

/•oo r /"OO p 

/ ada / dV = / bdb / d 4 k-- - , 

Jo J Jo Jc A or C v 

where on the right we choose any vector q £ C A with q 2 = b, and the /c-integral goes over C v 
or C A , depending on whether we started with a vector k on the lower or upper mass cone. 
Combining the contributions of the two summands in the square bracket in (|4.14p . we 
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obtain a fc-integral over C. Finally, the support property in (|4.6|) allows us to extend the 
fc-integration to all of M. We thus obtain 

5S = J bdbj Tr (A(k) P(k + q) JP(g)) + 5F(m 3 ,m 5 ) . 

Next, it is convenient to transform the /c-integral to a convolution integral, 

/ Wf A{k) P{k + q) = I M(k)P(q-k) =: (M*P)(q), 

where we flipped the sign of k and used the symmetry property of M in (|4.6|) . Then we 
can write the variation of the action in the more compact form 

5S = J Tr ((M * P)(q) 5P(q)^bdb + <5F(m 3 ,m 5 ) . (4.16) 

Our last step is to treat the variation of the function F. According to the chain rule, 

5F(m 3 ,m 5 ) = D 1 F(m 3 ,m 5 ) 5m 3 + D 2 F(m 3 ,m 5 ) 5m 5 . (4.17) 

It is an important observation that the relations (|4,4p determine M only up to singu- 
lar contributions on the light cone. More specifically, there is the natural freedom to 
modify M (£) by 5'- and ^-contributions of the form 

m(o ^ c £5'(f)e(e) + c^siexe). (4.18) 

As we shall see later (see Lemma 15.91) . by modifying Co and c± appropriately, we can 
precisely absorb the term (|4.17p which arises in the variation of F. Conversely, the freedom 
to choose the function F in (12. 12ft can be regarded as taking into account the fact that 
the distribution M is determined only modulo the contribution (|4. 18j) . In what follows, 
we shall always omit the term ^(rr^ms) in (I4.16|) . giving the following result. 

Proposition 4.2 The first variation of the action S, \2.12\) . can be written in momentum 
space as 

5S = 2 Tr (Q(q)SP(q)^bdb, (4.19) 
where q G C A is any vector in the lower mass cone with q 2 = b and 

Q(q) = \(M*P)(q) :=\j ^ M(p) P(q - p) . (4.20) 

Here M is the Fourier transform of the distribution M as introduced by h4-4\ )> f or specific 
values of the free parameters cq and c\ in \4.I8 ). 

To fully understand (|4.19l I4.2(J|) . one should observe that, according to (|2.ip . the fac- 
tor 5P(q) in (|4,19p is supported in C A . But for q in the lower mass cone, the convolution 
integral in (|4.20p is finite, because the integrand is compactly supported (see Figure 
The Fourier transform of Q is simply the product 

Q(0 = \m(0p(0- (4-21) 

The distribution Q(y — x) is precisely the kernel of the operator Q which appears in the 
EL equations [P, Q] = corresponding to the variational principle introduced in [H §3.5]. 

The notion of state stability imposes conditions on the eigenvalues of the matrix Q{q). 
For the sake of self-consistency we here restate the definition as given in [H Def. 5.6.2]. 
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p° / 


/ Pi 


/m(j>) 


*v' P(q-p)\ 

g 



Figure 4: The convolution Ai * P. 



Definition 4.3 The fermionic projector of the vacuum is called state stable if the cor- 
responding operator Q(k) is well-defined in the lower mass cone C A and can be written 
as 

U 

(4.22) 



Q(k) = a-^ + b 



(where \k\ 
properties: 



k 2 ) with continuous real functions a and b on C A having the following 

(i) a and b are Lorentz invariant, 

a = a {k 2 ) , b = b(k 2 ) . 

(ii) a is non-negative. 

(iii) The function a + b is minimal on the mass shells, 

(o + 6)(m|) = mf(a + b)(q 2 ) Va G {1, . . . ,g}. 



(4.23) 



In our setting, the representation (14.220 as well as (i) are obvious, and thus we need to 
analyze (ii) and (iii). The next theorem shows that the EL equations corresponding to 
the Lorentz invariant variational principle of Definition 12.21 give a necessary condition for 
state stability. Our proof will also explain the detailed form of the constraint (|2. 13|) . 

Theorem 4.4 The fermionic projector is a critical point of the variational principle of 
Definition if and only if the functions a and b in the representation ^4.22$ have the 
following properties: 

(a) (a + 6)K) = ( a + 6)(m|) Va, (3 G {1, . . . , g} 

(b) (a + 6)'(m«)= VaG{l,...,<7}. 
Proof. Using the ansatz (|2.ip . we find 

9 



m/3 



(4.24) 
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where q, qp £ C A are any vectors with q 2 = b and q^ = mk. Using (|4.22|) . this simplifies to 
1 r°° f \ 9 

- / Tr (Q(q) P{q)j b db = 2 £ pp mj (a + b)(m 2 ) . (4.25) 

Jo p=\ 

According to Proposition 14.21 5S is obtained from (|4,24|) by varying pp and mp, keeping 
the operator Q fixed. Equivalently, we can vary (|4.25p for fixed functions a and b. We treat 
the constraint (|2.13p as in (|2.17p with a Lagrange multiplier A. Varying the weights pp 
gives the conditions 

2mp{a + b)(m 2 p) = \mp V/3 € {1, . . . ,g} , (4.26) 
proving (a). Similarly, varying the masses rap and using (|4.26p gives (b). ■ 

The properties (a) and (b) in this theorem are weaker than the condition (iii) in the 
definition of state stability, because (|4.23p implies in addition that the points m/| are 
absolute minima of the function a + b. Moreover, our theorem makes no statement on 
the condition (ii). Nevertheless, the above theorem shows that our Lorentz invariant 
variational principle is closely related to state stability, and it seems a promising strategy 
for satisfying state stability to construct minimizers of our variational principle, and then 
to select those minimizers which satisfy the additional conditions of Definition 14.31 

In order to get stronger conditions from our variational principle, one needs to consider 
more general variations of the fermionic projector. A convenient method is to consider a 
"test Dirac sea" of infinitesimal weight and of possibly negative mass. Thus we consider 
instead of (|2.ip the fermionic projector 

9+1 

Hk) = J>/? (# + mp) 5(k 2 -m 2 p) Q(-k°) 

13=1 

and assume that p g +i vanishes for our unperturbed fermionic projector. We vary p g +i for 
given m g+ \ S M and introduce the function V by 

FK+l) = ^ dp^ Slp ^=° ■ 

We refer to V as the variation density. We point out that the {g + l) st Dirac sea should 
not be considered as being of physical significance; it is merely a mathematical tool to 
generate variations of the action. 

Theorem 4.5 The functions a and b in h4-2$fy are related to the variation density by 

V(m) = e(m) a(m 2 ) + b{m 2 ) . 
Proof. Taking into account that m g+ i may be negative, (|4.25p is modified to 
1 f°° f \ 9 

- / Tr (Q(q) P(q)j bdb = 2^9?™% {e(mp) a(m 2 ) + b{m 2 )) . 
Jo p=i 

Differentiating with respect to p g gives the result. ■ 

In view of this theorem, the conditions (ii) and (iii) are can be expressed in terms of the 
variation density by 

(ii') V(m) > V(-m) Mm G R+ 

(iii') V(mp) < inf R+ U V/3G {!,... g}. 
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5 Convolutions of Lorentz Invariant Distributions 



We saw in the previous section that it is useful to rewrite products in position space as 
convolutions in momentum space, 

FG(q) = (F * GO := J ^ F(k) G(q - k) , (5.1) 

and to analyze the resulting convolution integrals (see (|4.9p and (14.20P ). In this section 
we study such convolution integrals systematically. The distributions of interest like (13. 4p 
or (|3.13p are Lorentz invariant and are supported in the mass cone C. It is helpful to 
decompose these distributions into the sum of a distribution supported in the upper mass 
cone and a distribution supported in the lower mass cone. This gives rise to the following 
definition. 

Definition 5.1 A Lorentz-invariant distribution F G S'(M) is called negative if it is 

supported in C A and can be written as 

F(k) = f(k 2 )@(k 2 )Q(-k°) (5.2) 

with a real-valued distribution f G 5'(1R + ). Similarly, F(k) is called positive if F(—k) is 
negative. 

We remark for clarity that the representation (|5.2p poses an extra condition which in 
particular excludes all distributions supported at the origin. For example, the distri- 
bution 5 4 (k) is Lorentz invariant, supported in C A , but it cannot be represented in the 
form (pT2j) . 

Since F is completely characterized by /, it is most convenient to work only with the 
small letters. Note that for notational convenience, the function / does not carry a hat; to 
clarify momentum space we shall write out the arguments a := k 2 , . . .. In our computations 
we will use a compact notation which reminds of the corresponding operations in position 
space. The following list shows what operations on / and g mean for the corresponding 
distributions F and G: 

(f-_g)(a) (F*G)(k) ) 

/(a) ^ F(-k) \ (5.3) 

@f(a) i$F{k) ) 

The notation / cannot lead to confusion because /(a) is real-valued, and therefore the only 
reasonable meaning is to take the complex conjugate in position space, which corresponds 
to flipping the sign of k. If / is negative, then clearly / is positive. Keeping this in 
mind, we shall always assume that / and g correspond to negative distributions, whereas 
positive distributions are denote by / and ~g. The operator (j) is convenient for handling 
Dirac matrices. 

If F and G are both negative (or both positive), the integrand of the convolution 
integral (|5.ip is compact (see Figure [5] (left)). If on the other hand one factor is positive 
and the other is negative, the support of the integrand will in general be non-compact (see 
Figure [5] (right)). These two cases are essentially different and will be treated separately 
in Subsections 15.11 and 15.21 In Subsection 15.31 we shall analyze the distribution M, which 
already appeared in (|4.4l 14. 5p . 
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Figure 5: Convolution of two negative distributions (left) and a mixed convolution (right). 



5.1 Convolutions of Negative Distributions 

We introduce the function A(a, b, c) by 

A = a 2 + b 2 + c 2 - 2(ab + ac + be) . 
Note that A is symmetric in its three arguments. 



(5.4) 



Lemma 5.2 Suppose that f and g are negative distributions. Then the following convo- 
lutions are also negative and are given explicitly by 



(f'9)(a 
a(a 
(3(a 

{d k f-d k g){a 



1 



32vr 3 



dcf(c) 



dbg{b)^— 



1 



327T 3 







1 



32vr 3 
1 

32^3 



dcf(c) 

dcf(c) 
dcf(c) 







dbg{b)^ a -^p 

a + b — c 

~ c + b-a 
2a 



db 5 (6) VA 



o 



db g(b) V A 



(5.5) 

(5.6) 

(5.7) 
(5.8) 



Proof. Clearly, the convolutions are again Lorentz invariant, and a geometric argument 
(see Figure [5] (left)) immediately yields that they are again negative. To calculate the 
convolutions, we choose any q G C A with q 2 = a. Then 



(/•<?)(«) 



d*k 



(2vr)4 

JO 

del dbf(c)g(b) 



F(k) G(q-k) 



d 4 k 
(2vr) 4 



5(k 2 - c) Qi-k") S((q - k) 2 - b) e(k° - q°) 



poo poo p J4L 

= J dej^ dbf(c)g(b) J —6(k 2 -c) @(-k°)5(q 2 - 2qk + c-b)@(k°-q ). 

Due to Lorentz symmetry, we can assume that q points in the direction of k°. Introducing 
polar coordinates (u = k°,p = \k\, O G S 2 ), we can carry out the momentum integrals, 

-i poo poo 

(f-9)(a) = ^3 J dc J dbf{c)g{b) 
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fO 

X 



pV poo 

duj p 2 dp 5(oj 2 — p 2 — c) 5{a + 2ujyfa + c — b) 

Jqn JO 



la 

OO POO 



= 4tt3 Jo dc J db ^ c "> 9 ^ 9 ( u; o -c)X[-Va,o](^o) , (5.9) 
where x is the characteristic function, and u>q and po are determined by the conditions 

c + a — b 22 
U ° = -2y/E ' Po = W ° ~ ° ' ( ' 5 ' 10 ' ) 

A short calculation shows that j»q = A/ (4a) with A according to (|5.4p . Using that the 
Heaviside and the characteristic function in (15. 9ft vanish unless \/a > -y/c + \/&, we ob- 
tain (j53|) . 

In order to derive (15. 6h . we first compute 



oo poo p j4l 

dc jf d& /(c) g(b) J — vf i# 5{k 2 - c) 5(q 2 -2qk + c-b). 

We again choose a reference frame with q = (—y/a, 0) and introduce polar coordinates (u> = 
k°,p = \k\,ne S 2 ). Then 

-I POO POO 

m ■ g){a) = J dcj^ dbf(c)g(b)iu ^G(u 2 -c)l { _^ )0] (u )^ 

with ujq and po according to (|5.10p . Note that, due to spherical symmetry, the contribution 
~ 7 vanishes. A straightforward computation yields 

1 r a (<^-V~c) 2 n r/ n /iw ^ nx a-b + c v 7 ^ 
'o •/ 



(tf •*)(«) = =5-s / dc dbf(c)g(b) (-^7°) — ^ — = (W(a) 

327r J ./n In 2a a 



with a as in (|5.6p . The relation (|5.7p immediately follows from (|5.6p using the symmetry 
in / and g. 

To derive (|5.8p . we again assume that g points into the direction of k . Then, exploiting 
the spherical symmetry, 

(d k f ■ d k g)(a) = jS^ik F(k) i(q l - k l ) G(q - k) 

roo p Mh. 

dcj^ dbf(c)g(b) J ^—^(-kq + c)5{k 2 -a)5{q 2 -2qk + c-b) 

Y poo poo 



dc / db /(c) g(b) (uj y/a + c) G(u; - c) X[-y/E,o] ( w o) 



and a short computation gives (|5. 
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5.2 Mixed Convolutions 

We shall now analyze the convolution / ■ g of a positive and a negative distribution. We 
begin with the case q £ C A where the momentum is in the lower mass cone. This case is 
unproblematic because the integrand of the convolution integral (|5.ip becomes compact if 
we assume that f(k 2 ) vanishes identically for sufficiently large k 2 . Making this additional 
assumption, we can use the same method as in Lemma 15.21 



Lemma 5.3 Suppose that f and g are negative distributions and that f(k 2 ) vanishes 
identically for large k 2 . Then for q E C A and setting a = q 2 > 0, the following convolutions 
are well-defined and given explicitly by 

= ^-3/ dcf{c) dbg(b)^ (5.11) 



W-g)(a) = ($>a)(a), 

«(a) = ^sj a dc f^J dbg(b)VA (5.12) 
(/•$?) (a) = 

/?(«) = 3^3 ^ dc f( c )J o dbg(b)VAZ2—- (5.13) 

(d k f-d k g)(a) = dcf(c) dbg (b)VA^ . (5.14) 

32vr d J a J 2a 



(f'9)(a) = I 7K -^F(k)d(k-q) (5.15) 
dc I db /(c) g(b) I 5{k 2 - c) 0(-k°) S((k - qf - b) Q(q° - k°) 



Proof. The momentum vector q € C A satisfies the inequality q° < and thus 
d 4 k 

(27T) 

= jf dcjf dbf(c)g(b) J —5(k 2 -c)5(c-2qk + q 2 -b)e(q -k ). 

Due to Lorentz invariance, we can assume that q = (—1/0,6). Again choosing polar 
coordinates (uj,p, ft), we obtain 

-1 roo roo 

U-9){a) = dc J db /(c) 5 (6) 

du / p 2 dp 5(uj 2 -p 2 - c) 5(c + 2uy/a + a - b) 
-00 V 



1 

8^ 3 7o 



00 
00 



dc/ dbf{c)g(b) 

-\fa 



du Q(uj 2 - c) \Juj 2 - c <5(u; 2 - p 2 - c) <5(c + 2^^ + a - b) 

-00 



-00 
roo 

8vr 3 / ""'■'v-y^v-y y-'O " 



db /(c) 5(6) sJul-c—= 0(lu 2 - c ) 9(c - 6 - a) 



where we set 



a — 6 + c 

w ° = ^TT^" • 

-2^0 
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A straightforward calculation gives (|5.1ip . 

The relations (|5.12p - (|5.14p are obtained with the same method as the corresponding 
relations (|5.6p - (|5.8p in Lemma IQ1 ■ 

This lemma can also be applied in the case when q in the upper mass cone C v . Namely, 
changing the variables k — q — > q in (|5.15|) . one sees that flipping the sign of q corresponds 
to exchanging / and g. Hence the corresponding formulas are obtained if in (|5.11l 15. 14j) 
we exchange / and g on the right. 

We come to the most interesting case that q lies outside the mass cone, q 2 < 0. Since 
the hyperbolas {k 2 = c, k° < 0} and {(k - q) 2 = d , (k° - q°) < 0} (for c, d > 0) have a 
non-compact intersection, in this case the convolution integral (|5.ip necessarily diverges 
for any non-trivial distributions F and G. Our method to bypass this problem is to first 
break the Lorentz invariance by a regularization. This makes the convolution integral (|5.ip 
finite. Then we will subtract suitable counter terms which are supported on the light cone. 
After that, we can remove the regularization. We thus obtain Lorentz invariant formulas 
for the convolution integrals, which have to be understood "modulo infinite contributions 
supported on the light cone." This method of disregarding contributions on the light cone 
is motivated from the fact that we computed the pointwise product in (|2.3p only away 
from the light cone, with the intention of then extending the resulting function in the 
distributional sense across the light cone (|4.4p . The resulting distributions depend only 
on (|2.3p away from the light cone. Therefore, it is irrelevant to us how the regularized 
quantities behave on the light cone as the regularization is removed. This is why we can 
simply subtract divergent terms which are supported on the light cone. 

To regularize F, we set 

(f)(k) = f(k 2 )e £k ° (5.16) 
and introduce the distribution F £ S'(M) in analogy to (|5.2p by 

F £ (k) = f £ {k)Q(k 2 )@(-k°). (5.17) 

In this way, we can use the previous notation, with the only difference that the functions f e 
and g e are no longer Lorentz invariant and are thus no longer functions of q 2 , but of q 6 M. 
For simplicity, we assume that / and g have compact support. 

Lemma 5.4 Suppose that f(k 2 ) and g(k 2 ) are negative distributions which vanish identi- 
cally for large k 2 . Then for q 6 M \ C and setting a = q 2 < 0, the following formulas hold 
for the products of the corresponding regularized distributions \5. 1 6[ \5.17\ ), 

-i poo roc 

(f-¥)(q) = ^ dcf(c) dbg(b)H e (q,b,c) (5.18) 

1 f°° f°° h -4- r — n 
(d k f-d k W)(q) = dcf(c) dbg{b)—- H £ (q,b,c), (5.19) 

where the function H £ is given by 

H £ (q, b, c) = -j^ exp \e\q\ ^ + eq° . (5.20) 

Proof. Introducing the vector u = (£, 0), we can write the product f e g £ as 

f'OO roo 

(f-¥)(q) = / dcf(c) / dbg(b)I(q,b,d) 
Jo Jo 
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with 

b,c) = J 5(k 2 - c) 5((k - qf - b) 0(-k°) e(q° - k°) e ku+ ^ u . 

We choose a reference frame where q = (0, x, 0, 0) and u = (a, (3, 0, 0) with x > and a > 
Choosing cylindrical coordinates (u,p,r,ip) with k° = u> and k = (p, r cos ip, r sin tp), 
we can compute the integrals, 

^ f'Q roo roo 

I(q, b ,c) = ^IJ"IJPJ rdr 

x 6 (uj 2 -p 2 -r 2 -c) 8(u> 2 - (p - x) 2 -r 2 -b) e 2 — ^-*) 

D 2au-0{2p-x) 



r0 roo 

— duj dp Q(uj 2 -p 2 - c) 5(2px - x 2 + c - b) e 2 

«/ —CO «/ — CO 

- / dp 5{2px -x 2 + c-b) e -2«VP 2 + C -/3(2 P - a; ) 

^ ./ —CO 



16vr 3 
1 

32vr 3 
1 

647r 3 ax 



e -2aVK 2 +c-f3(2K-x) 



where we set 



x 2 — c + b 

K = 



2x 

This can also be written as 

I(q,b,c) = ^3 — exp (-ax A - fix B) , (5.21) 
where A and B are the expressions 

v/F+c V(-a- c + 6 ) 2 - 4oc v 7 ^ 

J~x — A 

x a a 

2K - x c-b 
B = = , 

x a 

which are covariant and independent of the regularization. Next, we transform back to 
the original reference frame where q = (q , q) and u = (e, 0). Since j3x = —uq = —eq° and 
ax = \J —u 2 q 2 + (uq) 2 = e\q\, we have the simple transformation rules 

ax -> e\q\ , (3x -> -eg . (5.22) 

Applying these transformation rules in (I5.21h gives (15.201) . 

The relation (|5,19p follows immediately from the transformations 



roo roo 

(d k f ■ d k f){q) = dc dbf(c)g(b)K(q,b,c) with 

Jo Jo 



d 4 k 

d 4 k 1 
(X) 1 2 

x<5(A; 2 - c) 5((jfe - g) 2 - b) Q(k°) e ku+ ^ u 
■- (a-c-b) I(q,b, c) . 



^ ((A; — (A; — g)) 2 - k 2 — (k — q) 2 ) 
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Obviously, the function H £ , (|5.20p . becomes singular as e \ 0. The next lemma shows 
that after dropping suitable terms which are supported on the light cone, the limit e \ 
exists and is Lorentz invariant. 

Lemma 5.5 Suppose that f(k 2 ) and g(k 2 ) are negative distributions which vanish identi- 
cally for large k 2 . Then the products of the corresponding regularized distributions $5.16\) 
have the decomposition 



(f e ■¥)(<!) 



Sl(q)+Rl(q) 
Sl(q)+Rl(q) 



where the Sf are distributions which are supported on the light cone, 

suppSf C {£ 2 = 0}, 
and the Rf are regular as e \ 0. The limits 

Ri 



lim R s : 



are the Lorentz invariant distributions 



Ri(q) 



dcf(c) dbg(b)Ki(q,b,c) 



with 



Ki(q,b, c) 



K 2 (q,b, c 



32vr 3 



<d{Vb 



Vb) 



\b-c\ 
a 

\b-c\ 



VE-\b-c\ 
2a 



1 



x < 



32vr 3 
y/A(b + c-a) 
2a 

VA(b + c-a) 



G(Vb 



2a 



e{^-^-Vb) 



\b- c| (b + c) 
2a 

\b-c\ (b + c) 
2a 



y/A (b + c-a)-\b-c\(b + c) 
4a 



6(6 -c) ifq£C v 
Q{c-b) ifqeC A 
ifq^C 

9(6 -c) ifqeC v 

9(c-6) ifqeC A 

tfq^C, 



where we again set a 



(5.23) 



(5.24) 



(5.25) 



It is interesting to compare the result of this lemma in the case q E C A with Lemma 15.31 
Then the first summands in (|5.24p and (|5.25j) coincide precisely with (|5.1ip and (|5.14p . 
respectively. The appearance of the additional second summands in (|5.24p and (|5.25|) is 
somewhat surprising; it shows that extending the convolution integral in the distributional 
sense to q € M yields additional contributions inside the mass cone. These additional 
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contributions can be understood from the fact they remove the poles of the convolution 
integral on the cone {q \ q 2 = 0}. Namely, using the expansion 

y/K = \b-c\+0(a) , 

one sees that K\ and K2 have no poles at a = 0. 

Proof of Lemma 15.51 We first consider the product f £ ■ g e . Extending the function H £ , 
(15,201) . by zero to the interior mass cone q £ C, we obtain the expansion 

H ^ c » - w + 4 nr 1 (c " b) + ^ e( - o) + 0(E) ' (5 ' 26) 

Here the poles on the mass cone have a mathematical meaning as distributional derivatives 
of the regular distribution 0(— a) log a. We implicitly use this procedure throughout the 
proof. The first two summands in (|5,26p are not Lorentz invariant, and we must analyze 
their behavior in position space. In preparation, we derive a useful formula for spherically 
symmetric distributions which are supported on the light cone. In position space, such 
a distribution can clearly be written in the form <?(£°) S(£, 2 )- We compute its Fourier 
transform in polar coordinates x = (t,r, tp), 



/(?) == / d A ig(e)Kf)e- iqx 

/oo roc rl 

dtg(t)e- iq0t / r 2 dr5(t 2 -r 2 ) / dcostfe^ 
-00 Jo J-i 



r cos ft 



2tt 



CO 



^ / r iq0t I rdr 5{t 2 - r 2 ) ( e^> - e~ l 

7T 

m 

vl / ,. ... ... _;„0-t / -L 



dtg(t) e-^ ( e *l«ll*l - e -*MI*l) 



dt g(t) e(t) (S^ - e"^') = — (%° + \q\) - h(q° - \q\)) , 

Ao\ J-oc V J 2|g| 

where /i is the Fourier transform of the function 2mg(t)e(t). Since the function g can 
be arbitrary, the function h is also arbitrary. Also, the above calculation clearly carries 
through if h is a one-dimensional distribution. We conclude that for any distribution h, 
the expression 

m = ^(% + k1)-%°-k1)) (5.27) 

is supported on the light cone. Choosing h(x) = e(x), one sees that the distribution 
2Q(—a)/\q\, and thus also the first summand in (|5.26p . are supported on the light cone. 
The second summand is not of the required form. But choosing h(x) = (c — b) e(x)/(4x), 
one sees that its extension 



J w 2a 



is supported on the light cone. Setting 



-l ifgeC v 

1 if q G C A 



©(-«) \ \b~c\ PP 
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S £ is supported on the light cone, and the difference H e — S e has a Lorentz invariant limit, 



lim(tf e - S £ 

e\0 



\b-c\ 
a 

\b-c\ 
a 

VA- \b-c 



2a 



9(6 -c) if q£C v 
@(c-b) ifqeC A 



If q £ C A , we must also take into account the contribution from Lemma 15.31 In the case 
q € C v , we can apply Lemma [5731 after double conjugation, 

dcs(c) / dbf{b)^—. 



(f-g)(Q) = (9 ■/)(-?) 



32vr d 



Adding all these contributions gives the result for the product f £ ■ g e . 

In order to treat the product c^/ 6 • <9 fc <? £ , we note that (|5.14|) and (|5.19p differ from 
the corresponding formulas (|5.1ip and (|5.18|) only by the additional factor (6 + c — a)/2. 
Since b and c were treated as fixed parameters throughout the proof, this additional factor 
can be written as the differential operator in position space (b + c + □?)/2. Using that 
the distributions D^Sf are again supported on the light cone, the above arguments all go 
through after applying the operator (6+c+D^)/2. From this consideration we obtain (|5.23p 
with 

K 2 (q,b,c) = Ki(q,b,c) . 



Finally, we subtract the contribution 

r \b-c\ 



2 

16 — c| 



9(6 -c) ifgeC v 
9(c-6) ifqeC A 



(5.28) 



\b-c\ 



liq^C 



which for any fixed 6 and c can be written as c\ + C2 @(q 2 ) e(<?°) with constants c\ and 
The Fourier transform of the first summand is a multiple of the distribution 5 4 (£). Using 
that \3l(G(q 2 ) e(q )) = in the distributional sense, one sees that the second summand is 
supported on the light cone (alternatively, see Corollary 16.21 below) . We conclude that the 
contribution (15.28H is supported on the light cone {£ | £ 2 = 0}, completing the proof. ■ 

This lemma shows that the products / -g and d^f ■ d k g are singular only on the light cone. 
Furthermore, it gives a useful explicit formula for the regular contribution away from the 
light cone in momentum space. Note that the regularization is needed in order to make the 
singular contribution on the light cone finite, but the remaining regular contribution does 
not depend on the regularization. With this in mind, we write the result of Lemma 15.51 in 
the compact form 

(f-g)(q) = (l.c.)+/ dcf(c) dbg^K^qAc) (5.29) 
Jo Jo 

f'OO f'OQ 

(d k f ■ d k g)(q) = 0-C-)+/ dcf(c) db g{b) K 2 (q,b,c) , (5.30) 

Jo Jo 
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where "(I.e.)" stands for a singular (possibly infinite) contribution on the light cone. 
The next lemma treats the products which involve one derivative. 

Lemma 5.6 Suppose that f and g are negative, Lorentz invariant distributions. Then, 
using the short notation introduced before $5.29\) . 



W-9)(q) = (l.c)+/ dcf(c) dbg{b)Li{q,b,c) 

Jo Jo 

f'OO f'OQ 

(f-m(d) = (l-c)+/ dcf(c) dbg(b)L 2 (q,b,c) 

Jo Jo 



(5.31) 
(5.32) 



with 



Li(q,b, c) 



x < 



32vr 3 

a — b + c 

2a 2 
a — b + c 



2a 2 



Q(Vb- y/E- v^)+ {b f a2 2ab Q(b-c) i/?£C v 
e(V~c -^-Vb)- {h ~ f~ 2ab 9(c - b) z/<7GCM5.33) 



2a 2 



. a — b + c 

A + 

4a 2 



2ab 



4a 2 



ifq$C 



L2(q,b, c) 
( 



x < 



32vr 3 

a + b — c 

2a 2 
a + b — c 

2a 2 



Q(Vb - Va" - Vc~) - {b C n )2 9 2aC 6(b-c) i/geC v 



e(v^-v 

a + 6 — c 
4a 2 



v / 6) + 
i — c 



2a 2 

(6 - c) 2 - 2ac 



2a 2 



2 -2ac 



4a 2 



6(c-6) z/gGC A (5-34) 



where again a = q 2 . 

Proof. Using that differentiation in position space corresponds to multiplication in mo- 
mentum space and vice versa, we obtain from (15 . 16H that 



2\ „efc u 



\f £ {k) = 4f(k 2 )e 



y k (h(k 2 )e^° 



ie 7 °) h £ (k) 



where the negative, regularized distribution h £ is defined by 



h £ {k) = h(k 2 )e ek ° 



and h{k 2 ) = / /(e) de . 
Jo 



Now we can compute the product <f)f ■ g with the help of Lemma 15.5 



1 



$f-g = lim-(-£-isry )h £ gz 



(l.c.)+/ dch(c) I dbg(b)-±K x (q,b,c) 
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Applying the transformations $ = — ift and 

roc roc rc poo roc 

/ dch(c) Ki(q,b,c) = dc de f(e) K(q,b,c) = / de /(e) / dcKi(q,b,c), 
Jo Jo Jo Jo Je 

we obtain the representation (|5.3ip with 

Li(q,b,e) = J -ft q Ki{q,b,c) dc . 

In order to compute L\ more explicitly, we substitute for K\ the expression (|5.24p 
and carry out the c-integration. The following indefinite integrals are easily verified by 
differentiation, 



/ 



d VA 
da a 
d \b-c\ 
da a 

The asymptotic expansion for large c 



dc 
dc 



b — c — a 



\b-c\ 



2a 2 
b-c 



2a 2 



v/A = (c-6) - a^±| + 0(c~ 2 ) 
c — o 

^(b-c-a) = -(c-bf + 2ab + a 2 + 0{c~ l ) 
shows that the difference of the two indefinite integrals has a finite limit as c — > oo, 



. b — c — a ., . b — c 
hm | VA — — \b-c\ 



2a 2 



2a 2 



b 1 
a + 2- 



Taking into account the boundary conditions as determined by the Heaviside functions 
in (|5,24p gives for L\ the expression 



32vr J 



2a 2 
a — b + c 



2a 2 



S(Vc- Va~- Vb) 



2a 2 
(b-c) 2 



0(6 -c) 



if q € C x 



a — b + c |6 — c| (6 — c) b 1 
4a 2 + i^ 2 + 2^ + 4 



if q & C . 



We finally subtract the distribution 



PP 1 1 

bQ (b-c)— + -Q(a)e(-q°) + -, 

which can be written in the form (|5.27p and is therefore localized on the light cone. This 
gives (|5.33p . 

The formula for L2, (|5.34p . is obtained from (|5.33j) by double conjugation, f-ftg = ftg-f- 
This gives the identity L,2(q, b, c) = Li(—q, c,b). ■ 

It is instructive to verify the Leibniz rule 

Kf-9)-W-g)-(f-m = 
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by substituting the results of Lemma 15.51 and Lemma 15.61 This gives 



with 



/•oo 1*00 

(Kf ■ 9) - W -9)-(f- ffi) (?) = (l.c) + / dc /(c) / db g(b) M(q, b, c) 

Jo Jo 



M(q, b, c) = ty Ki (q, b, c) - L x (q, b, c) - L 2 (q, b, c) 

t^e(6-c) + ^e(6- c ) \i q ec y 



32vr J 



a 

\b-c\ 
a 



2a 

-G(c-6)- — e(c-b) if qeC A 
a 

\b — c\ b-c , .„ , n 



5.3 Convolutions Involving Dirac Seas 

We now apply the previous results on convolution integrals to composite expressions in 
the fermionic projector. We begin with the convolution integral (|4.10p . 

Theorem 5.7 For any k G C with k 2 > 0, 



i t 



64 vr 3 k 4 



e ( k °) ^2 PaP/3 J(k 2 ,m a ,m/3) 
a,p=l 



(5.35) 



where 



J(a, x, y) = -v/A(a,x 2 ,y 2 ) (x - y) e(x 2 - y 2 ) [{x + yf - a] @((\x\ - \y\) 2 - a) (5.36) 
with A according to J5.^| ). 
Proof. Introducing the distributions 



/(a) = ^pp5{a 
(3=1 



ffl< 



and g(a) = ^ppmp 5{a - m|) , (5.37) 

/9=1 



we can write the fermionic projector (|2.1|) with the notation (|5.3p as P = —iftf + g. Then 



(p-p)o = -i((m-9-g-(m)- 



(5.38) 



For & G C A , the appearing convolution integrals were computed in Lemma 15.31 We thus 
obtain 

(F-P)o = -i?a(a) 

x (e(v^- >/o- v^)-G(v / 6 
where we set a = k 2 . Using (|5.37p . we can carry out the integrals to obtain (|5.35p with 

J(a,x,y) = -2x V / A(a, x 2 , y 2 ) (a - x 2 + y 2 ) (0(|y| - ^ - \x\) - 6(|x| - y/a - \y\)) 
= -2x ^&{a,x 2 ,y 2 ) (a - x 2 + y 2 ) e{\y\ - \x\) 9((|y| - |x|) 2 - a) . 
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As (|5.35p is symmetric in a and (5, we may symmetrize J in x and y, giving (|5.36|) . 
For q G C v , the result follows similarly using the remark after the proof of Lemma 15.31 ■ 

Next we can introduce the distribution Ai, which already appeared in our sketch as the 
Fourier transform of the distribution Ai defined by (|4.4p . Here we proceed in the opposite 
way: we first define Ai and then show that its Fourier transform really satisfies (|4.4p . We 
denote this Fourier transform without a tilde by Ai. This clarifies that Ai is a specific 
distribution satisfying (|4,4p , whereas Ai is not unique but can be changed freely according 

to flngp . 

Theorem 5.8 The function 

1 U 9 

M(k) = ^3 ^0(fc 2 )e(fc°) £ p a Pp{J(k\m a ,m p ) + K(k 2 ,m a ,m p )) (5.39) 

a,/3=l 

with J according to \5. 36\) and 

K(a,x,y) = (x-y)\x + yf - 2a(x 3 + y 3 ) (5.40) 

defines a tempered distribution, if the pole on the cone {k \ k 2 = 0} is understood as the 
distributional derivative of a logarithm, i.e. 

^e{k 2 )e(k°) stands for - @ k (log(k 2 ) Q(k 2 ) e{k )) . (5.41) 

Its Fourier transform Ai(£) satisfies away from the light cone the relation 

M(0 = 2A (0 y^L (5.42) 

with Aq as given by \2. 3\) . 

Proof. It suffices to consider any pair a, (3 G {1, . . . ,g}. In the case m a ^ mp, we write 
out (|5.42p in momentum space to obtain similar to (15.38P the condition 

M{k) = -2t((0/). g-g-m) ■ 

Computing the convolution integrals with the help of Lemma 15.61 modulo singular con- 
tributions on the light cone, a straightforward calculation similar to that in the proof of 
Theorem 14.41 gives (|5.39p with K according to (|5.40p . Expanding the factor \fK in (|5.36p 
for small a, 

y/A(a,b, c) = \b-c\-a + 0{a 2 ) , (5.43) 

one sees that Ai{k) has no poles on the cone {k \ k 2 = 0}. Thus Ai(k) is a L^-function, 
which clearly defines a tempered distribution. 

In the case m a = mp, the function J vanishes identically, and thus 

M{k) = -±^^ € (k°)e(k 2 )p a ppml. 

Interpreting the pole according to (|5.4ip . this is again a tempered distribution. It has 
such a simple form that (|5.42p can be verified by explicit calculation. Namely, the Fourier 
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transform of Xi can be given in closed form (see Corollary 16.41 below). Comparing this 
formula with an explicit computation of Aq in position shows that (|5.42p is indeed satisfied. 



Let us briefly discuss the result of Theorem 15.71 and Theorem 15.81 We first point out that, 
due to the Heaviside function in (|5.36j) . the function JF{k) vanishes identically for large k 2 , 

T(k) = if k 2 > max^gj! . j9 } . (5.44) 
However, it has poles as k 2 \ 0, as one sees from the expansion (|5.43p . 



* $e(fc°) ^ ( {m a -mp) 2 {m a + m p f mj + m 3 - 



a, (3 with m a ^m,fj 

+}^O{k ). (5.45) 

The distribution A4 differs from T by the function K in (15, 39ft and by the fact that the 
poles of M have a meaning in the distributional sense (|5.4ip . The order of the pole is 
smaller than that in (]5.45p . 

M(k) = -^#&(k 2 )e(k ) P-PP \2 ~ + ¥O(k ). (5.46) 

a,P with m a =ma 

Since (|5.35p vanishes identically for large k 2 , using (|2.9l 12. 10p we find that M. has the 
simple asymptotics 

M(k) = 2tt 2 £ e(k 2 ) e(k°) + ^ if k 2 > max^ {1 ,..., 9} m 2 . (5.47) 

We end this section by specifying the effect of the parameters c$ and c\ in (|4.18p on 
the integral in (|4.16p . 

Lemma 5.9 The Fourier transforms of the distributions 

A(0 := t5(e)e(e), B(0 := ^'(£ 2 ) e^ ) 

satisfy the relations 

Tr ((i* P){q) 6 P(q)^ a da = 32vr 4 <5m 3 (5.48) 
J Ti^(B*P){q)5P(q)^ada = -32^ 4 5m 5 , (5.49) 

where rri3 and rri5 are defined by \2.9\ \2.10\) . and where again a = q 2 . 

Proof. Again using the notation (|5.3p . we write the fermionic projector as P = —i<j)f + g 
with / and g according to (I5.37p . The Fourier transforms of A and B can be computed 
explicitly (see Corollary I6.2I below). Decomposing them into a sum of a positive and a 
negative distribution, we write A and B as 

i$(h(a) - h(a)) 
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with 

h(a) = 8iv 2 5'(a) fori and h(a) = 2vr 2 5(a) for B . 
The convolutions P * A and P * B are then expressed by 

+ ig) ■ {pi - 



For any q G C A , we can write these convolutions with the help of Lemma 15.21 and 
Lemma \52 



($)f-{ph-ph)){a) = a(a), 

a{a) = ^Jo 

(ig-(0h-0h))(a) = i#9(o), 

/3(a) 



327T 3 Jo 



dc f(c) e(a — c) 



dc g(c) e(a — c) 



db h(b) VA 



dbh(b) VA 



c + b — a 
2a 

a + b — c 
2a? 



We now compute the remaining integrals using the special form of the functions /, g 
and h. We begin with the calculation for A. Using similar to (I5.43P the expansion 



VA = \a-c\-b 



a + c 



a — c 



+ 0(b 2 



(5.50) 



we obtain 



a(a) = — — / dc f(c) e(a — c)—i V A 



1 
1 

47T 



r/ 



db 
1 6 



a=l 



(5(a) 



1 

47T 

1 

47T 



dc g(c) e(a — c 



dcg(c) 



o 



r/6 

I b 



c + 6 — a 
2^ 



a + 6 — c 



6=0 



2a 2 



6=0 



5 3 

4^^- Pa a 2 

a=l 



and thus 



(P*A)(q) = "^E ^ W + a) • 
o=l 

Now we take the trace with a "test Dirac sea" 

R(q) := p(i + m)8(q 2 -m 2 ). 

This gives 

f°° / \ 1 9 

/ Tr ((P*i)(g)%)) ada = Vpp a K + m 3 ). (5.51) 

V 1 ^fi 

After varying p and m, we can build up <5P by taking the sum of g test Dirac seas. Using 
the symmetry of (|5.5ip in p a and ttt- , fn a ^ wc thus obtain 

r-&((P*4)( 9 )5P(g))ado = ~si^p a p p (m 3 a + m 3 ^, 
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and comparing with (|2.9p gives (|5.48|) . 

The calculation for B is similar. Using (|5.50p . we obtain 

a(a) = 7^-/ dc f '(c) e(a - c)y/A ° 



2a 1 6=0 



\2 1 _JL ^2 _ ^2 

167T /n 2a 327r 



a=l 

2 1 9 ™ ^2 „\2 



1 f 00 (c-a) 2 1 ^ m a (m£ - a) 



16vr Jo yw 2a 2 32vr 

a=l 



and thus 



(p *£)(<?) = -^-ffe K n2 fl)2 M + a ) 

OZ7T tt 
ct=l 

J Tr \(P* B)(q) R(q)J a da = -— ^ pp Q (m Q - m) 2 (m a + m)' 



a=l 



Comparing with (|2.10p gives (j5.49j) . 



6 A Lorentz Invariant Regular izat ion and its Fourier Trans- 
form 

According to (|4.4p and (|2,15l I2.8P , on the light cone M. has the leading pole 

M(0 = m 3 le(e)e(C ) + 0(r 2 ), (6.1) 

and this pole needs to be regularized in order to give the action a mathematical meaning. 
In this section we shall construct a Lorentz invariant regularization of the pole in (|6.ip . 
which is explicit both in position and in momentum space. Our analysis is based on the 
following Fourier integral. 

Lemma 6.1 The following identity holds in the sense of distributions: 

I 6( " 2) e(k0) ^ = " W ^ ^ + 2^ 6{e) £ «°> (6 ' 2) 

Proof. Due to spherical symmetry, we can assume that £ = (t, r, 0, 0) with r > 0. Choosing 
polar coordinates k = (w,pcos$,psin?9cos(^,psin$sin(^), the Fourier integral becomes 

A := f ^ e - eJ t9(k 2 )e(k )e ik t 



(2vr) 4 

1 /*°° rl^l / /*i 

/ aW(u/)e^* / p 2 dp e~ ^ ~ p2 ) / dcostf e" ipr cos1? 

due^e** pdpe-^ 2 - p2 ) (e~ ipr - e ipr ) 



8ir 3 r 
i 

87r 3 r 



00 



due(u) pdpe-i^ ~ p l+^t-pr) 
-00 J— \u\ 
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We now introduce the "mass cone coordinates' 

1 



2 (w+p) , v = - (u-p), 



and light cone coordinates 



s = - (t-r), I = - (t + r) 



Then our integrals transform to 



47T 3 ? 



JO 



oo <j — oo 



du dv (u — v) e 



Carrying out the v-integral in the case u > gives 



(u — v) e 



-2euv+2ivl 



dv 



o 



i d 
U+ 2di 



-2euv+2ivl 



dv 



- 2euv+2ius+2ivl 



i d 







1 



2 <9Zy 2en - 2iZ 

whereas in the case it < we get the same expression with the opposite sign. Hence 

1 



A 



47T 3 ? 



i d 



2 dl 2eu - 2il 



e 2ius du 



\ ( d d 



87r 3 r V ds dl 



d d 



i 



2eu - 2il 



e 2ius du 



87r 2 re \ ds dl 



e(s) e(s/)e— r 



where in the last step we computed the integral with residues. We finally transform back 
to polar coordinates, 



A 



1 d 



e{t) e(t 2 - r 2 ) e — *r 



4tt 2 e r dr 

carry out the distributional derivatives and write the result in Lorentz invariant form. ■ 

Before going on, we collect of a few simple Fourier integrals which can easily be obtained 
from (IQl. 



(6.3) 
(6.4) 
(6.5) 
(6.6) 



Corollary 6.2 The following equations hold in the sense of distributions: 

£'(£ 2 M£°) e~* fc? A = -i7r 2 Q(k 2 )e(k°) 
5{£ 2 )e(f) e- ik Zd 4 £ = -Am 2 5(k 2 ) e(fc°) 
J $6'{£ 2 )e(£ ) e~ ik S d 4 £ = 2ir 2 k 1 5(k 2 ) e(k°) 
J £<5(£ 2 )e(£°) e- ,fc$ A = 8vr 2 ^^(A; 2 ) e(A; ) 
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Proof. The identity (|6.3p could be obtained by a direct computation of the Fourier inte- 
grals. We here use another method where we recover (|6.3p as a limiting case of Lemma [6. 11 
Using that 

2 „ 1 



1 

£ Jo 



~J2 ' e 2e dz 



e JO 



z e 2 £ 



one sees that 



J_ e -^e(z) - lj(z) ^ 



2e 



-4<5'(z) (6.7) 
with convergence in <S'(R). This also shows that (|6.2p converges in S'{M) to 



/ 



(27T)' 



e(fc 2 ) e(fc°) e' 



On .ifc^ 



7T 



5'(e 2 )^°), 



except maybe at the origin ^ = 0, where the light cone has its cusp. To prove convergence 
at the origin, one verifies easily by a scaling argument that for any test function r/ £ 

C§°(M), 



lim lim / r? 

<5\o e\o ' 



ie-ie(£ 2 )e(£°) 



2s 



. 



(6.8) 



We finally apply Plancherel's theorem to (|6.7p to obtain (|6.3p . 

Applying the operator to (|6.3p gives (|6.4p . To derive (|6.5p . we apply the opera- 
tor i(f) k to (|6.3|) . Finally, the identity (|6.6|) follows from (|6.4|) by applying the differential 
operator i(j) k . ■ 



The main importance of Lemma 16.11 lies in the fact that integrating over e gives the 
Fourier transformations of regularized poles. 

Proposition 6.3 The following equations hold in the sense of distributions: 



d A k 

WT 4 



^n k (log(k 2 )Q(k 2 )e(k°))e ik Z 



1 

d^k 



7T 



9(e 2 )e(e°) + <5(£ 2 )e(£°)(c-l + loge) 



(2tt)' 



^0 k (log(k 2 )e(k 2 )e(k )) 



0\\ 



? (l- e -fe]eK 2 )e(f°) + Hf) e(«°)(c + log e) 



Here the constant c is given by 



and 7 is Euler's constant. 



c = 7 - log 2 , 



(6.9) 

(6.10) 
(6.11) 



Proof. Integrating the variable e in (|6.2p over the compact interval [e, L] with L > e gives 
the identity 



/• d A k 2 

y (^o 1 ^ 



e 2 



e"^- e(fc 2 )e(fc°)e^ 



i 2 



4vr 2 £ 2 



T7T e 2_L 



e-fe )e(C 2 )e(e°) + ^<5(£ 2 )e(£ )(logL-log £ ) 
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Bringing the term involving logL to the left, we can use (|6.4|) to combine it with the 
integrand, 



d A k f 1 



(2tt) 4 [k 2 



ek 2 Lk' 1 

e 2 — e 2 



9(fc 2 )e(fc°) - 5{k 2 ) e(k°) (c + logLUe 



i 

"top 



1 

e 



A direct calculation yields for any r\ 6 Cq°(1R) and x > the relations 
inn (^j r](a) ^ ((1 - e - ^) 0(a)) da - logL^) 



(6.12) 



?7(0) (log(s) + 7 - log 2) + 



17(a) 



1 d 



a da \ da 
Thus, choosing c according to (|6.1ip . we obtain 

lim < - (e~^ — e~^\ 9(a) — 5(a)(c — 1 + logL) 



a^ (9(a) log a) ) = r/(0) (log(a;) + 1) + 



(6.13) 
(6.14) 



-- W(log(a) 0(a)) e _ T- 



where W is the differential operator (|3.15p . and the derivatives and the convergence are 
meant in <S'(M). Using at the origin a scaling argument similar to (|6.8f) . we can thus take 
the limit L^ooon the left of ([6TT2]) to obtain (pJUj) . 

In order to derive (16.10p . we apply the operator — ifi^ to (|6.12p . Using (I6.13|) together 
with the identity 



v( a ) ~r (©( a ) l°g a ) = ^(O)log(a;) + o(x), 
da 



we conclude that 



lim 

L— >oo I a 



{!( e_¥ " e ~ ¥ ) 0(a) " 5(a)(c + logL) } e(A;O) = lh{^g(a)e(a)e(k°))e-^ , 



giving (j6.10p . ■ 

The method which we just used to take the limit L — > 00 can be also be applied to take 
the limit e \ in (|6.9p and (|6.10p . giving the following result. 

Corollary 6.4 The following equations hold in the sense of distributions: 

d 4 k 



(2vr) 4 
d A k 



a k (log(k 2 )Q(k 2 )e(k )) 



0\\ Jk£ 



@ k (log(k 2 ) Q(k 2 )e(k )) e ifc « 



1 

tor 2 
1 

1^2 



— n^io g (e)e(e)e(e)) (6.15) 

^((log(e 2 )-l)9(^ 2 )6(e )) (6.16) 



Let us explain why the previous Proposition is so useful for regularizing Ai near the 
light cone. 9ur goal is to introduce a regularization M £ which is explicit both in position 
and momentum space and makes our action finite, i.e. 



{M £ ) 2 ada < 



(6.17) 
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(see also (|4.1|) and the discussion thereafter). Carrying out the derivatives, one sees that 
the distribution (|6.16p has a pole of the desired form (|6.ip . 

Thus (|6.16p extends this pole in the distributional sense across the light cone and gives an 
explicit formula for the corresponding Fourier transform. But this does clearly not solve 
the problem of making the integral (|6.17j) finite. The integrand in (|6.10j) differs from that 

in (|6.16p by the factor e~~ . The decay of this factor for large k 2 can be thought of as a 
"smooth momentum cutoff," and thus we refer to it as the 

regularizing factor e~~ in momentum space. (6.18) 

The interesting point is the effect of this regularization on the pole in position space as 
described by the right side of (|6.10p . The term <5(£ 2 ) leads to a distributional contribution 
supported on the light cone. It is not obvious how to make sense of it in an integral of the 
form (|6,17p . But we can simply drop this contribution by considering instead of (|6.17p an 
integral of the form 

poo roc 

/ (M £ ) 2 ada := lim / (M £ ) 2 ada (6.19) 

(we remark that an equivalent, but maybe a bit cleaner method for getting rid of the 
distribution supported at a = is to subtract this distribution from (|6.17p and to use that, 
according to (16. 5|) . this changes the Fourier transform only by an irrelevant contribution 
supported on the mass cone {k \ k 2 = 0}). Disregarding the ^-distribution, the effect of 
the regularization can be described by the 



replacement 
Using the formula 



J_ _d_ 

z 2 dz 



1 



in position space. (6.20) 



1 / * \ f 2e 



^ l-e-Zj = J e ~ Xz dX 

one sees that the square bracket in (|6.20p is a smooth function in the variable z. Also, 
due to the decay of the factor e~2?, the regularization is indeed localized in a strip z ~ e 
around the light cone. With (|6.18p and (|6.20p we have a simple regularization method, 
which is convenient both in momentum and position space. 

We are now ready to complete the constructions outlined in Section HJ According 
to the left of (|6.16p . the pole ~ corresponds in momentum space to a term of the 
form lfc/k 2 . Thus the relevant term for the regularization of A4 can be read off from the 
expansion (|5.47p . 

Definition 6.5 For any e > and k G C with k 2 > 0, we introduce the function J\f £ by 
M £ (k) = -Ll e (/c°) p a ppK{k 2 ,m a , mf5 ) + 2^ 2 m 3 ^ (e~^ - l) e(k°) , 

n a, 13=1 ^ ' 

where K is the function |5.^0[ ). For any e > 0, we define the distribution M £ by 



I 

k' z 

with M. according to \5. 39\) . The Fourier transform of M. £ is denoted by M. £ . 



M £ {k) = M(k) + 2ir 2 m 3 -j^ [ e~^ - l) S(k 2 ) e{k°) (6.21) 
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Proof of Lemma \4-l\ From the expansion (|5.46p we know that for any e > 0, the 
pole of Ai £ (k) on the mass cone is integrable in the L 2 (M, a da)-norm. Furthermore, 
from (|5.47p and Definition 16.51 we see that M £ (k) for large k 2 has the form 

M £ (k) = 2vr 2 ©(A; 2 ) e(k°) e'^ + if fc 2 > max^^ s} m 2 . (6.22) 

We conclude that 7W £ G L 2 (M,ada), and thus Corollary 13.21 applies. 

jf Tr((X £ ) 2 )zdz = M(M £ ) 2 )ada, 

were we disregard the distributional contributions at z = and a = as explained 
before (|6.19p . From (|6.22|) we can also determine the counter term needed in (|4.7|) . Namely, 
substituting (|6.22p into (|4.7p and computing the integral, we find that choosing 

F e (m 3 ,m 5 ) = G(m 3 ,m 5 ) - 4?r ^ m 3 _ 32vr 4 m 3 m 5 log e , 

the limit e \ exist. The limit clearly coincides with that in (|2.12p for a suitable choice 
of the function F — G. 

The identity (|4.9p follows immediately from Definition 16.51 an d thus it remains to jus- 
tify the limits in (|4.1ip and (|4.12p . Since T vanishes for large k 2 according to (|5.44p . the 
integral clearly converges for large a. For the behavior near a = 0, we use the expan- 
sion (|5.45p to rewrite the integral in (|4.1ip as 



c + 0(a))da, (6.23) 



Im £ 

a 

where c stands for the leading coefficient in (|5.45p . According to Definition 16.51 and (|5.39|) . 
the square bracket converges in 5'(M), and therefore (|6.23p makes sense in the limit e \ 
if we interpret the round brackets as the test function. 

To analyze (|4.12p . we split up the integral into an integral over the interval a G (0, 1] 
plus an integral over [1, oo). Using the explicit form of M s , the integral over (0,1] can 
again be written in the form (|6.23p . which establishes convergence as e \ 0. The inte- 
gral over [1, oo), on the other hand, can be analyzed exactly as that in (|4.7j) . because Af £ 
and M £ coincide for large a according to (|4.9l 15.44)) . ■ 



7 Additional Free Parameters 

In view of the application to state stability, the following consideration gives rise to ad- 
ditional free parameters. As is shown in [U Appendix A], it is possible to regularize 
the fermionic projector in such a way that the property of a distribution .M-P-product 
is no longer satisfied, but where the continuum limit yields an additional distributional 
contributions supported at £ = 0. More precisely, (|4.2ip is replaced by 

Q(0 = \M(Z)P(0 + c 2 5\0-c z ^5\0-caU5^) 
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with real parameters c 2 , C3 and C4, which can be chosen arbitrarily. Let us consider how 
these additional parameters can be built into the variational principle of Definition 12.11 
Taking the Fourier transform, we see that the additional terms change the functions a 
and b in the representation (|4.22p by 

a(k 2 ) -» a{k 2 ) + c 3 \k\ , b{k 2 ) -» b(k 2 ) + c 2 + ci k 2 . (7.1) 

The parameter c 2 changes Q only by an irrelevant constant, and thus we may disregard 
this parameter. Then, going through the calculation in the proof of Theorem 14.4^ one 
immediately verifies that the terms in (|7.ip can be arranged by modifying the action as 
follows. 

Definition 7.1 We define the extended action S ex t by 

9 9 
Sext = S + c 3 y^^ppm 4 + c 4 y] J o 3 m| 
0=1 0=1 

with S as in \2. 6\) . Here C3 and C4 are two free real parameters. The corresponding 
variational principle with constraint \2.13\i is referred to as the extended variational 
principle. 



8 Numerical Construction of Minimizers 

We first rewrite the action in a form which is most convenient for the numerical analysis. 

Proposition 8.1 (Extended action in momentum space) We introduce the real- 
valued function 

H(a,x,y) = - (j(a,x,y) + K(a,x,y)) (8.1) 



a, 

(with the functions J and K as given by 15.361 \5.40j )) and choose a parameter 

a max > max ml . 
/3e{i,..,s} p 

Then the extended action (see Definitions \7jj\ and \2.1\) can be written as 

Sext = 216^10 ^ PaPpP^Ps / H(a,m a ,mp) H(a,m^,m 8 ) da (8.2) 

9 9 

+ F(a max ,m 3 ,m 5 ) + c 3 ^ pp m% + c^^ppm^. (8.3) 

0=1 0=1 



Proof. It clearly suffices to consider S, because the extended action is obtained simply by 
adc 
M( 
by 



adding the extra terms in Definition 17. 11 From (|6.2ip and (|5.47p it is clear that A4 £ (k) £ —>° 
M(k) uniformly on the set {k\ < k 2 < a max }, whereas for large k 2 , M e is given explicitly 



M e (k) = 2n 2 f Q(k 2 ) e(k°) e~^ + ^) if k 2 > a max 
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Figure 6: The variation density in the case of one generation of mass m\ = 1 and the 
parameters cq = C3 = C4 = and c\ = 2 • 10 6 (left), c\ = 10 7 (right). 



Hence, using the result of Theorem 15.81 in the formula for the regularized action in mo- 
mentum space (|4.7p . we obtain 



1 ramax 

<S = 16 1Q ^ PaPpP-yPS / H(a,m a ,mp) H(a,m y ,m s ) da 

+ J5^e\fl ( 47f4 / ( m 3 e ~^~ +4m 5 a -1 ^ da + F e (m 3 ,m 5 )^ . 

The result follows because the last line depends only on a max , ITI3 and m.5. ■ 

Note that the argument before (|5.46j) shows that the function H, (18, lj) . is bounded near a = 
0. Hence in (|8.2p we only need to integrate a bounded continuous function over a compact 
interval. This makes it possible to compute the action with standard numerical methods. 
Using standard tools of nonlinear optimization, one can construct numerical minimizers. 

Alternatively, one can construct critical points by searching for numerical solutions 
of the conditions (a) and (b) in Theorem 14.41 It is helpful that the resulting equations 
are polynomials in the parameters pg, making it possible to use specific algorithms for 
computing the roots of systems of polynomials. For a detailed description of the numerics 
we refer to the PhD thesis [8]; here we proceed by discussing a few of the numerical 
solutions. 

In the case g = 1 of one generation, we constructed all minimizers of the unextended 
action (i.e. C3 = c = 4 = 0). By scaling we can assume that mi = 1. We found a one- 
parameter family of minimizers parametrized by ci, whereas cq must vanish. In Figure [6] 
the variation density is shown in typical examples. In agreement with Theorem 14.41 and 
Theorem 14. 5[ we see that V'(l) = 0. However, the solutions are not state stable, because 
the condition (iii) in Definition 14.31 (or equivalently the condition (hi') on page fT7|) is 
violated. More specifically, the variation density has near m = 1 the expansion 

V(m) = (m- l) 2 (aci +^log|m- 1|) + o((m - l) 2 ) 

with constants a, (3 > 0. Therefore, it is impossible to arrange that m = 1 is a local 
minimum of V. Since the extra contributions of the extended action are polynomials 
in m, this statement remains true for the extended action. We conclude that there is no 
state-stable Dirac sea configuration with one generation. In words, our action principle 
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prefers to "split up" the Dirac sea into two seas. However, by choosing c\ sufficiently large 
we can arrange that the minima of V are arbitrarily close to m = 1 (see Figure [6] right). 

The case g = 2 of two generations gives some insight into the nature of our variational 
principle. By scaling we can always assume that m\ = 1 < 7772 and p\ = 1. Thus for 
the unextended action we have the four free parameters 777,2, P2, Co an d c\ to satisfy the 
three conditions V'(m\) = = V (777,2) and V(mi) = ^(7772). Therefore, one expects 
for given 7772 a discrete number of minima. Indeed, we found for every 7772 exactly one 
minimizer. In Figure \7\ the variation density of the minimizers is shown for different values 
of 7772- If 7772 is close to mi, the result can be understood in analogy to the minimizers 
for one generation: again the points m = rrij are local maxima of the variation density. 
Thus the solution is not state stable, and our action principle has the tendency to split 
up the two Dirac seas into even more seas. However, if 777,2 is increased, a "potential wall" 
emerges between m\ and 7712, such that the number of local minima of V increases from 
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three to four (see the cases m<i = 4.25 and mi = 4.5). The points m = rrij are still local 
maxima, but the neighboring local minima are very close (see the case mi = 5). If mi is 
further increased, the local minima approach the points m = rrij exponentially fast, and 
thus the resulting variation density looks very much like a state stable configuration (see 
the case mi = 10). Keeping in mind that the physical masses of the three generations 
of leptons scale with factors of the order 100, and that a relative uncertainty in the rest 
masses of 10 -100 can certainly not be measured in experiments, in such cases it is fair to 
identify the points rrij with the neighboring minima, so that we may regard the rrij as 
absolute minima of V. 

In the case g = 3 of three generations, we can assume by scaling that 1 = mi < 
m 2 < m 3- We found different classes of minimizers. In view of the physical masses of the 
leptons, we were most interested in constructing solutions where the masses are far apart, 
i.e. m\ <C m2 <C m%. Our attempts to construct such solutions were successful only for 
the extended action. Figure [8] gives an example of a state stable Dirac sea configuration 
with three generations. The corresponding values of the parameters are mi = 5, m% = 20, 
Pl = 1, p 2 = 10" 4 , p 3 = 9.696 • 10" 6 , c = -6.692 • 10 8 , a = -2.516 • 10 9 , c 3 = 9658.25, 
c 4 = 8416.56. 

As already mentioned in the introduction, the systematic study of the minimizers goes 
beyond the scope of this paper. Clarifying the detailed structure of the minimizers and 
critical points of our action is an important project for the future. Talking of future 
developments, we finally point out that information on the constants Co, ci, ci and C4 
could be obtained by studying regularized or discretized Dirac sea systems (for example 
using numerical methods [7] for large systems). Moreover, the weight factors pp have an 
influence on the continuum limit (see Appendix A]), and thus the detailed analysis of 
the continuum limit for interacting systems should give constraints for the weight factors. 
In view of these future investigations, which should reduce our number of free parameters, 
the introduced action principle seems to be a promising step towards a physical theory 
which makes predictions for the ratios of the masses of elementary particles. 
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